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Abstract. The yield of muon pairs in the invariant mass region 1<M<2.5 GeV/c 2 produced in heavy-ion 
collisions significantly exceeds the sum of the two expected contributions, Drell-Yan dimuons and muon 
pairs from the decays of D meson pairs. These sources properly account for the dimuons produced in 
proton-nucleus collisions. In this paper, we show that dimuons are also produced in excess in 158 A GeV 
In-In collisions. We furthermore observe, by tagging the dimuon vertices, that this excess is not due to 
enhanced D meson production, but made of prompt muon pairs, as expected from a source of thermal 
dimuons specific to high-energy nucleus-nucleus collisions. The yield of this excess increases significantly 
from peripheral to central collisions, both with respect to the Drell-Yan yield and to the number of nucleons 
participating in the collisions. Furthermore, the transverse mass distributions of the excess dimuons are 
well described by an exponential function, with inverse slope values around 190 MeV. The values are 
independent of mass and significantly lower than those found at masses below 1 GeV/c 2 , rising there up 
to 250 MeV due to radial flow. This suggests the emission source of thermal dimuons above 1 GeV/c 2 to 
be of largely partonic origin, when radial flow has not yet built up. 



PACS. 14.40.Lb - 25.75.Nq - 25.75.Cj 

1 Introduction 

The intermediate mass region (IMR) of the dimuon mass 
spectrum, between the <j> and the J/ip resonances, is ex- 
pected to be well suited to search for thermal dimuon 
production [1], due to the favourable relative production 
yield with respect to the other contributions (Drell-Yan 
dimuons, meson decays, etc). 

a Corresponding author: ruben.shahoyan@cern.ch 



Intermediate mass dimuon production in proton-nucle- 
us and heavy-ion collisions was previously investigated by 
NA38 [2] and HELIOS-3 [3] in p-W and S-U(W) collisions 
at 200 GeV, and by NA50 [4, 5] in p-A (where A stands 
for Al, Cu, Ag, W) and Pb-Pb collisions, respectively at 
450 and 158 GeV. All three experiments reported a very 
reasonable description of the IMR opposite-sign dimuon 
mass continuum measured in the "elementary" proton- 
nucleus collisions. This continuum could be accounted for 
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Fig. 1. Schematic representation of the NA60 experimental layout. 



by a superposition of the two expected "signal" processes, 
the Drell-Yan dimuons and the muon pairs resulting from 
simultaneous semi-muonic decays of (correlated) D and 
D mesons, on the top of a "background" contribution. 
The latter is due to uncorrelated decays of pions and 
kaons and can be estimated from the measured like-sign 
muon pairs. In contrast, all experiments observed that the 
opposite-sign dimuon samples collected in the heavy-ion 
collision systems, taking into account the estimated back- 
ground contribution, significantly exceeded the level ex- 
pected from Drell-Yan and "open charm" sources, calcu- 
lated using the same procedures that successfully repro- 
duced the p-A data. 

As discussed in [6] , two prime interpretations were able 
to describe the findings of NA38 and NA50 equally well: 
the long-sought thermal dimuons [7], and an increase in 
the open charm production cross section per nucleon, from 
p-A to A-A collisions [8]. Several other reasons which 
could increase the yield of IMR dimuons were also pro- 
posed. Alternatively to charm enhancement, the number 
of muon pairs entering the phase space window of the 
experiment could be increased, while the total charm pro- 
duction cross section remains unchanged. This could re- 
sult from, e.g., the smearing of the D/D pair correla- 
tion resulting from rescattering of the charmed quarks 
or mesons in the surrounding dense matter [9]. It was 
also suggested that the mass spectrum of the Drell-Yan 
dimuons could be modified in heavy-ion collisions with re- 
spect to the proton-nucleus case, because of higher-twist 
effects that increase the yield of low mass dimuons [10]. 
Another source of dimuons that could be present in heavy- 
ion collisions, while being negligible in p-nuclcus colli- 
sions, is secondary Drell-Yan production, where the quark- 
antiquark annihilation uses valence antiquarks from pro- 
duced pions [11]. 

A decisive step in understanding the origin of the ex- 
cess dimuons is to clarify the decade-long ambiguity be- 
tween prompt dimuons and off- vertex muon pairs. In the 
first case they can be thermal dimuons or extra Drell-Yan 
dimuons; in the second case they result from decays of D 
mesons, which have a relatively long lifetime: ct — 312 /im 
for the D + and 123 for the D°. The clarification of the 



physical origin of the IMR dimuon excess was one of the 
main motivations of the NA60 experiment. Thanks to its 
ability to measure the offset of the muons with respect to 
the interaction vertex, NA60 can separate, on a statistical 
basis, the prompt dimuons from the off-vertex muon pairs. 

In this paper we present a study of the intermediate 
mass dimuons produced in In-In collisions at 158 GeV/nuc- 
leon, based on data collected by the NA60 experiment 
in 2003. The paper is organized as follows: Section 2 de- 
scribes the NA60 experimental setup, the data reconstruc- 
tion procedure and the general performance of the appa- 
ratus; Section 3 explains in some detail the background 
subtraction procedure; Section 4 presents the results. Pre- 
liminary results were presented before [12]. 



2 The NA60 experimental setup and data 
reconstruction 

Figure 1 shows a general view of the NA60 apparatus. Its 
main components are the muon spectrometer (MS), previ- 
ously used by the NA38 and NA50 experiments [13], and a 
novel radiation-hard silicon pixel vertex tracker (VT) with 
high granularity and high readout speed [14] , placed inside 
a 2.5 T dipole magnet just downstream of the targets (see 
Fig. 2). The first spectrometer is separated from the sec- 



dipole, 2.5 T 



beam tracker 



I 



tot 

box vertex 

HHHI 

tracker 



10 cm 
i 1 



I 
I 

IC 



muon 
filter 



Fig. 2. View of the target region. From the left: 2 stations of 
the beam tracker followed by 7 Indium targets and 16 planes 
of the vertex tracker. 



ond by a hadron absorber with a total effective thickness 
of ~ 14 Aint and — 50 X . 
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Before entering into more details, we list the key fea- 
tures of this somewhat unique setup: 

— The vertex telescope tracks all charged particles up- 
stream of the hadron absorber and determines their 
momenta independently of the MS, free from multiple 
scattering effects and energy loss fluctuations in the 
absorber. The matching of the muon tracks before and 
after the absorber, both in coordinate and momentum 
space, strongly improves the dimuon mass resolution in 
the low-mass region (less so higher up), significantly 
reduces the combinatorial background due to n and 
K decays and makes it possible to measure the muon 
offset with respect to the primary interaction vertex. 

— The additional bend by the dipole field in the target 
region deflects muons with lower momenta into the 
acceptance of the MS, thereby strongly enhancing the 
opposite-sign dimuon acceptance, in particular at low 
masses and low transverse momenta, with respect to 
all previous dimuon experiments. A complete accep- 
tance map in two-dimensional M-p T space is contained 
in [15]. 

— The selective dimuon trigger and the radiation-hard 
vertex tracker with its high read-out speed allow the 
experiment to run at very high rates for extended peri- 
ods, maintaining the original high luminosity of dimuon 
experiments despite the addition of an open spectrom- 
eter. 

A detailed description of the muon spectrometer can 
be found in [13]. Its magnet defines the rapidity window 
where the dimuons are accepted, 3 < j/i a b < 4. The trigger 
system is based on four hodoscopes along the beam di- 
rection. Each of them has hexagonal symmetry, with the 
sextants operating independently The system provides a 
highly selective dimuon trigger requiring that the four ho- 
doscope slabs hit by each muon match one of the prede- 
fined patterns. This ensures that the muon was produced 
in the target region. In addition, the trigger imposes that 
the two muons must be detected in different sextants. 

A detailed description of the radiation-hard silicon pixel 
vertex tracker used by NA60 in 2003 can be found in [16]. 
The readout pixel chips, developed for the ALICE and 
LHC-B experiments [17], operate with a 10 MHz clock 
frequency. The VT can provide up to 12 space points per 
track, 9 of them from planes oriented such that the hori- 
zontal coordinate (in the bending plane) is measured with 
higher precision. 

The target system is composed of seven Indium sub- 
targets, 1.5 mm thick each and separated by an inter- 
distance of 8 mm, adding up to an interaction probabil- 
ity of 16% for the incident Indium ions. An interaction 
counter is located downstream of the VT. It is made of 
2 scintillator blades, appropriately holed to let the beam 
pass through. They are independently read by 2 photo- 
multiplier tubes in coincidence and thus allow to tag the 
interactions in the target region. A beam tracker [18], 
made of two tracking stations 20 cm apart, was placed 
upstream of the target system. The beam tracker allows 
to measure the flight path of the incoming ions and to de- 
rive the transverse coordinates of the interaction point, in 



the target, with an accuracy of 20 /xm, independently of 
the collision centrality. A zero degree calorimeter, previ- 
ously used in NA50, is located in the beam line, inside the 
muon filter, just upstream of the Uranium beam dump. It 
estimates the centrality of each nucleus-nucleus collision 
through the measurement of the energy deposited by the 
beam "spectator" nucleons. 

Around 230 million dimuon triggers were recorded on 
tape during the 2003 run. In approximately half of these 
events a dimuon was reconstructed from the muon spec- 
trometer data. Two data samples were collected, with dif- 
ferent currents in the ACM toroidal magnet: 4000 A and 
6500 A. The higher field reduces the acceptance in the 
highly populated region of low transverse mass (hit) muon 
pairs, thereby increasing the number of high mass dimuon 
events collected per day, for a constant lifetime of the data 
acquisition system. The details of the event selection can 
be found in [16]. 

The reconstruction of the raw data proceeds in several 
steps. First, muon tracks are determined from the data of 
the eight MWPCs and validated by the hits recorded in 
the trigger hodoscopes. If the event has two reconstructed 
muon tracks which fulfil the trigger conditions, the tracks 
in the silicon planes of the vertex tracker are also recon- 
structed, and the interaction vertices are searched for. The 
track reconstruction efficiency is ~ 95 % for peripheral In- 
In collisions and ~ 90 % for the most central ones. The key 
step in the data reconstruction is the matching between 
the muon track, extrapolated from the muon spectrome- 
ter to the target region, and the charged tracks found in 
the vertex tracker. This is done by selecting those associa- 
tions between the MS tracks and the VT tracks which give 
the smallest weighted squared distance (matching \ 2 ) be- 
tween these two tracks, in the space of angles and inverse 
momenta, taking into account their error matrices. The 
matching procedure combines the good MS momentum 
resolution (a p /p ~ 2%) with the excellent VT angular 
precision (~ 1 mrad) to obtain the kinematics of the muon 
before undergoing the multiple scattering and energy loss 
induced by the hadron absorber. This procedure, as men- 
tioned before, improves the dimuon mass resolution, from 
70-80 MeV/c 2 to 20-25 McV/c 2 in the u and <j> mass re- 
gion, and allows to correlate the muon's trajectory with 
the interaction vertex, the point of primary interest for 
the present paper. 

The resolution of the vertex determination, and its de- 
pendence on the number of tracks associated with the ver- 
tex, can be obtained from the dispersion between the mea- 
surements provided by the beam tracker and by the vertex 
tracker. As shown in Fig. 3, it is better than 10 ^m in x 
and 15 /xm in y, except for the most peripheral collisions. 

Figure 4 shows the distribution of the z coordinate of 
the reconstructed vertices, for the events with only one 
vertex found. We can clearly identify the seven In targets, 
placed between the two windows that keep the target box 
in vacuum, downstream of the two beam tracker stations, 
also placed in vacuum. 

The resolution of the offset distance measured between 
the matched muon tracks and the collision vertex can be 
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Fig. 3. Dispersion between the transverse coordinates of the 
interaction vertex given by backtracing the VT tracks and by 
extrapolating the beam tracker measurement (open symbols), 
as a function of the number of tracks attached to the vertex. 
Derived vertex resolution (solid symbols). 
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Fig. 4. Distribution of the z coordinate of the interaction ver- 
tex for the events that only have one good vertex reconstructed. 



evaluated by using the muons from 3 ftp decays. Since the 
B — > J/ip contribution is negligible at SPS energies, all the 
J/ip mesons are promptly produced. Moreover, the pion 
and kaon decay background is negligible under the J/t/j 
peak. Therefore, the offset distribution made with these 
muons (shown in Fig. 5) directly reflects the resolution of 
the muon offset measurement: 37 nm in x (bending plane) 
and 45 /um in y. These values are the convolution of the 
track uncertainties with the accuracy of the transverse 
coordinates of the vertex. 

Since the offset resolution of the matched tracks is 
affected by multiple scattering in the silicon planes, the 
analysis is performed using a weighted muon offset vari- 
able, essentially insensitive to the particle's momentum. 
Its definition is 
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Fig. 5. Offset distribution for matched muons from J/ip de- 
cays, in the x (circles) and y (triangles) coordinates, normalized 
to unit area. 

where V -1 is the inverse error matrix accounting for the 
uncertainties of the vertex fit and of the muon kinematics 
fit. Ax and Ay are the differences between the coordinates 
of the vertex and those of the extrapolated muon track, in 
the transverse plane crossing the beam axis at z = ^vertex- 
We then characterize the dimuon through the weighted 
dimuon offset, defined as 

A„„ -- 
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3 Background subtraction 



'AxW^ 1 



AyWyy 1 + 2AxAyV a 



xy 



(1) 



The sample of collected dimuons include a combinatorial 
background (CB) originating from decays of uncorrelated 
hadrons, mostly 7r's and K's. Since the production and de- 
cay of the parent hadrons are independent processes, this 
combinatorial background (CB) contributes in the same 
way to the opposite-sign and like-sign samples of muon 
pairs, increasing quadratically with the charged particle 
multiplicity. It is important to emphasise that the NA60 
trigger treats in exactly the same way the opposite-sign 
and the like-sign muon pairs. 

While this kind of background was already present in 
previous dimuon experiments, such as NA38 and NA50, 
a new type of background appears in the NA60 dimuon 
spectra due to the matching procedure. Indeed, any VT 
track having a small enough matching \ 2 with respect 
to the muon track is considered a matching candidate. 
For high enough charged particle multiplicities, the muon 
track will have several possible matches and, naturally, at 
most one of them is correct. All the other associations are 
fake matches and, if selected, need to be subtracted. Pairs 
where one muon match is fake and the other one is cor- 
rect are also part of the fake matches dimuon background 
(FB). Because the probability of having a fake match, for 
each muon, is proportional to the track density, the yield 
of matched muon pairs where only one muon is fake rises 
linearly with the charged particle multiplicity, while the 
yield of doubly fake pairs rises quadratically. 
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Fig. 6 shows the single muon matching \ 2 distributions 
for correct and fake matches, normalized to unit area. The 
distribution for fake matches is obtained from the real 
data using the mixed- events technique (see section 3.2): 
each muon from the MS is matched with the VT tracks 
of another event with the same vertex position and mul- 
tiplicity. Since the probability of having a fake match is 
determined only by the density of the non-muon tracks 
in the phase space of the matching parameters, the ob- 
tained distribution of the fake matches is automatically 
obtained with correct normalization. By subtracting it 
from all matches, the distribution for the correct ones is 
obtained. One can see that the distribution for the fake 
matches is much natter than the corresponding distribu- 
tion for the correct matches. Selecting exclusively matches 
with a matching \ 2 below a certain threshold value, the 
signal-to-background ratio can be improved at the expense 
of losing some fraction of the signal. 
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Fig. 6. Matching \ 2 distributions for correct and fake matches, 
normalized to unit area. 

Despite generating the fake matches background, the 
muon matching procedure leads to a very significant re- 
duction of the combinatorial background due to muons 
from pion and kaon decays. This happens for two rea- 
sons. First, when the pions or kaons decay within the ver- 
tex tracker, the kink at the decay point prevents most 
of the muon tracks from being reconstructed by the sil- 
icon tracker, thereby removing these decay muons from 
the matched muons sample. Second, if the pions or kaons 
decay downstream of the vertex tracker, the matching be- 
tween the meson track and the muon track usually results 
in rather large matching \ 2 values due to the kink be- 
tween the parent meson and decay muon. By only select- 
ing matches of relatively low \ 2 values, we suppress the 
yield of correctly matched muons from 7r,K decays. 



3.1 Combinatorial background 

The combinatorial background contribution to the oppo- 
site-sign dimuon distributions can be evaluated from the 



measured like-sign muon pair samples, using an event mix- 
ing technique. Two muons from different events with simi- 
lar characteristics are randomly picked and paired to build 
the "mixed CB" sample. Only muon pairs respecting the 
dimuon trigger conditions are kept (in particular, the mu- 
ons must be in different sextants) . Additionally, a sextant- 
dependent weight is applied to the measured muons to cor- 
rect the bias introduced by the "different-sextants" trigger 
requirement. The mixing is done separately in 40 narrow 
bins in ccntrality in order to avoid the bias due to the 
variation of the single muon kinematics and the /z+//Lt_ 
ratio within the bin width. The technical details arc given 
in Appendix A. It is important to notice that this mixed- 
event procedure reproduces not only the shape of the CB 
spectra but also their absolute normalization. 

An important issue to consider in the CB generation 
concerns the selection of the muons used for the mixing 
and for the computation of the sextant-dependent weights: 
they must share the same target, the same field polar- 
ities, the same charged track multiplicity bin, and the 
same configuration of the silicon pixel planes (same ac- 
ceptances, efficiencies, etc). Since we are interested in de- 
termining the background contributions to the matched 
opposite-sign dimuon spectra, it could seem natural to 
only use for the event mixing single muons which have at 
least one match. This would ensure that the normaliza- 
tions, given by Eq. (A-7), would be correctly computed 
for the matched dimuons spectra. Unfortunately, in the 
case of the analysis presented in this paper, the event 
mixing must be made using all muons (including non- 
matched muons), for reasons explained in the next sec- 
tion. Therefore, the computed normalizations correspond 
to all the dimuons reconstructed in the muon spectrome- 
ter, matched and non-matched. 

The accuracy of the mixed-event method can be evalu- 
ated by comparing the mixed and measured distributions 
of like-sign muon pairs, as shown in Figs. 7 and 8, respec- 
tively for the mass and weighted offset (see also section 3.4 
below) . 



3.2 Fake matches background 

The most direct way to evaluate the fake matches back- 
ground is to use the overlay Monte Carlo method: the VT 
hits of generated dimuons are superimposed on the VT 
hits of measured data, in order to ensure realistic occu- 
pancy conditions. By definition, every matched track built 
without the sufficient number of the Monte Carlo muon 
hits is a fake match (conventionally we define the match as 
fake if the track has more than two non-muon hits). This 
method works well for studies that do not use the muon 
offset information. However, while the kinematic variables 
of the dimuons are quite robust with respect to the un- 
avoidable differences between the measured and simulated 
data (in particular the residual misalignment of the geo- 
metrical setup), the offset distribution is much more sen- 
sitive. Therefore, for studies which require the muon offset 
information, such as the study reported in this paper, we 
had to develop a fake subtraction procedure which only 
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uses measured data and provides reliable spectra in any 
dimuon parameter. 

This is again achieved by an event mixing technique, 
matching the muons of a given event with the VT tracks 
of other events (with the same production target, charged 
multiplicity, running conditions, etc.), for the real data 
and for the artificial CB sample. Since the real data sample 
also contains dimuons where only one of the two muons 
was incorrectly matched, we added to our simulated FB 
sample a certain proportion of muon pairs where one of 
the muons retains its own match (see Appendix B). 

In principle, there could be two alternative approaches 
for reproducing the matched dimuons. If the two muon 
tracks from the MS have mi and m.2 matches, then there 
are mi x m 2 matched dimuons. Given the different shapes 
of the x 2 distributions for correct and fake matches (Fig. 6), 
the pair composed by the two best matches (those with 
the smallest \ 2 ) has the highest probability of being the 
correct one. However, this probability decreases with the 
increase of the number of matching candidates, leading 
to a degradation of the correct matching efficiency from 
peripheral to central collisions. Also, as explained in Ap- 
pendix B, since it appears impossible to estimate analyti- 
cally the amount of FB in such best matches spectra, one 
should rely on the overlay Monte Carlo method. Alterna- 
tively, one can consider all mi x m 2 matches for a given 
MS dimuon. Although this will increase the amount of FB 
to subtract - hence increasing the statistical error - such 



11 all matches dimuon spectra" offer a better control of the 
systematic errors, and allow us to cross-check the Monte 
Carlo based subtraction method. The results presented in 
this paper are obtained using this "all matches" proce- 
dure. Finally, to reproduce the offsets of the fake matches 
at the interaction vertex, we apply the same algorithm as 
used for the CB (see Appendix B). 

The top panel of Fig. 9 shows the mass distribution of 
Monte Carlo uj dimuons (with "all matches") and those 
of the fake contributions, both as estimated by Monte 
Carlo tagging and by event mixing. The bottom panel 
shows the ratio between the mixed and the MC-tagged 
fake backgrounds. Both methods are seen to agree in the 
mass range having meaningful statistics. 

3.3 Extraction of the correctly matched signal 

The CB and FB contributions obtained via event mixing 
are not independent from each other. As explained above, 
the combinatorial pairs can be correct or fake matches, 
CB = CB corr + CBf a kc. Clearly, the latter is also present 
in the total sample of "fake" pairs, due to signal or to 
combinatorial dimuons, FB = FB S i g + FB c b, with FB c b 
= CBf a kc. Therefore, if we subtract the total fake back- 
ground and the total combinatorial background from the 
data, we do not obtain the correctly matched signal, be- 
cause the "fake combinatorial" pairs are subtracted twice: 
Data - FB - CB = Data - FB si „ - CB corr - 2CB fakc = 
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Signal — CBfako- To avoid this double subtraction, we must 
estimate the "fake combinatorials" contribution. This is 
done by applying the CB estimation algorithm described 
in Appendix A to the generated FB sample, using its like- 
sign pairs for the event mixing. 

Figure 10 shows the mass spectra (integrated over all 
collision centralities) for the measured opposite-sign dimu- 
ons, full combinatorial and fake signal background sources, 
and the extracted correctly matched signal. 

Figure 11 shows the dimuon weighted offsets distribu- 
tions. It is worth noticing that the region 2 < < 6, de- 
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Fig. 10. The measured opposite-sign dimuon mass spectrum 
(points), the signal (bullets: correctly matched; solid line: in- 
correctly matched), and the combinatorial background (both 
correctly and incorrectly matched: dashed line) for the 4000 A 
data with matching x 2 <l-5. 



cisive to disentangle the open charm contribution from the 
prompt one, has a significant level of incorrectly matched 
signal. 



3.4 Systematic errors from the background subtraction 

The yield of like-sign dimuons remaining after the subtrac- 
tion of the mixed-event spectra constitutes a very good es- 
timate of the residual combinatorial background contribu- 
tion left in the opposite-sign spectra. Based on the ratios 
shown in Fig. 7, 8, the accuracy of the generated back- 
ground is estimated to be ~1%. The resulting system- 
atic uncertainty of the extracted signal is defined by the 
signal-to-background ratio which strongly depends on the 
kinematic bin and the cut imposed on the matching \ 2 - 
Being comparable to the statistical error, it changes from 
— 25% for p T <0.25 GeV/c at masses around 1.2 GeV/c 2 
to —1% near 2.5 GeV/c 2 with a lose cut x 2 =3, and it 
improves by more than a factor 2 for a x 2 =1.5 cut. 

To account for these systematic errors in the several 
fits mentioned in the next section, whenever necessary, the 
statistical errors of the estimated background were glob- 
ally scaled up to ensure that the residuals of the like-sign 
spectra are compatible with zero within three standard 
deviations. Even in the worst case (when fitting a pt dis- 
tribution in a narrow dimuon mass range), the errors were 
not increased by more than — 10%. 

The results presented in the next section were checked 
by repeating the fits with both cuts on the % 2 and both 
low and high magnetic field data and were found to be 
always compatible within the quoted errors. 



4 Results 

4.1 Expected sources of IMR dimuons 

The analysis was separately performed for the 4000 A and 
6500 A event samples. Dimuons were selected in the kine- 
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Fig. 11. Data, signal and background spectra for weighted 
dimuon offset distributions (defined in Eq. 2) in IMR for 4000 A 
data with matching \ 2 <1.5. 
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matical domain denned by < y cms < 1 and | cos 9cs I < 
0.5, where #cs is the Collins-Soper decay angle. Only events 
with a single reconstructed vertex in one of the seven In- 
dium targets were kept for the analysis. Since the signal- 
to-background ratio strongly depends on the matching \ 2 
cut (see Fig. 6), the data were analyzed with two different 
cuts (x 2 < 1.5 and \ 2 < 3) to evaluate possible systematic 
effects related to the background subtraction. 

The shapes of the Drell-Yan and open charm con- 
tributions were obtained with the Pythia Monte Carlo 
event generator [19] version 6.325, using the CTEQ6L 
set of parton distribution functions [20] including nuclear 
effects through the EKS98 model [21]. The primordial 
fcx was generated from a Gaussian distribution, of width 
0.8 GeV/c for Drell-Yan (to describe the pt distribution 
of dimuons heavier than the J/VO and 1.0 GeV/c for the 
open charm [22]. The charm quark mass was kept at the 
default value of 1.5 GeV/c 2 and the Py parameter (proba- 
bility that the c-quark hadronizcs in vector states) was set 
to 0.6 [23], resulting in an effective cc — > fi + fi~ branching 
ratio 0.84%. 

The Pythia events were generated at the vertices of 
real events and propagated through the experimental setup 
using the Geant 3.2 transport code [24]. The resulting hits 
were added to those of the real event used to set the in- 
teraction vertex and the resulting overlay Monte Carlo 
events were then reconstructed, with the codes used to 
process the measured data. We only kept the events sur- 
viving the selection cuts also applied to the real data. It 
is worth noticing that the | coshes I < 0.5 window signifi- 
cantly cuts the open charm contribution, because Pythia's 
strong D/D pair correlations give a coshes distribution 
peaked at —1 and +1. This means that the fraction of 
accepted dimuons from D pair decays is small and very 
sensitive to the kinematic distributions and correlations 
used in Pythia. Figure 12 shows the cos 9qs distribution 
of the cc — * dimuons, with mass in the range 1.16— 

2.56 GeV/c 2 , before and after the reconstruction step. 




Fig. 12. cos#cs distribution of cc — > fi + y,~ dimuons in the 
1.16-2.56 GeV/c 2 range, as generated by Pythia, before (upper 
line, normalized to unit area) and after (lower line) applying 
the muon reconstruction algorithm. 



The basic goal of the analysis is to compare the mea- 
sured yield of signal IMR dimuons to the expected yield, 
from the Drell-Yan and open charm contributions. Since 
this comparison is to be done as a function of collision cen- 
trality, the events are distributed in 12 sub-samples, using 
the number of tracks reconstructed in the VT as centrality 
estimator. The corresponding average number of nucleons 
participating in the In-In collision, N pai t, can be derived 
within the Wounded Nucleon Model N c h ~ qN paTt , where 
_/V c h is the number of charged particles produced in the 
VT angular window, corrected for acceptances and effi- 
ciencies [25]. The proportionality constant q is found by 
matching the observed dN/dN^ to dN / dN palt generated 
from the Glauber model. 

To judge if the sum of the expected contributions re- 
produces or not, in amplitude and shape, the measured 
mass and weighted offset IMR dimuon signal distribu- 
tions, we must determine the normalizations of the Drell- 
Yan and open charm spectra. We start by fixing the rel- 
ative normalization of the open charm contribution with 
respect to the Drell-Yan contribution, by calculating the 
ratio of their production cross sections. In order to cal- 
culate this ratio, we use the high-mass Drell-Yan cross 
sections measured by NA3 [26] and by NA50 [27], which 
show that Pythia's Drell-Yan spectrum needs to be up- 
scaled by a K-factor of 1.9, and we use a charm cross 
section of a cS — 8.6 ph. The latter number corresponds to 
the value required to reproduce the dimuon mass distribu- 
tions measured by NA50 in p-A collisions at 450 GeV [4, 
5], 36 ±3.5 fib, scaled down in energy (to 158 GeV) using 
Pythia. It is worth noting that this value is around 1.8 
times higher than what would be obtained from a "world 
average" estimate [22] based on measurements of D meson 
production from fully reconstructed hadronic decays. One 
should stress that neither NA50 nor NA60 are suited to 
measure the full phase space cc cross section: as can be 
seen from Fig. 12, the acceptance window | cos6*cg| < 0.5 
defined by the muon spectrometer contains less than 20 % 
of all dimuons from DD decays, and the extrapolation of 
the measurement to full phase space strongly depends on 
how well the correlations between the two decay muons are 
described by Pythia. The obtained open charm to Drell- 
Yan cross-section ratio is then kept the same for all cen- 
trality bins, since both processes are expected to scale with 
centrality in exactly the same way (proportionally to the 
number of binary nucleon- nucleon collisions) . 

The next step is to determine the Drell-Yan normaliza- 
tion from the yield of high-mass dimuons. Given the rela- 
tively small statistics, especially when the events are sub- 
divided in several centrality bins, we calculate the Drell- 
Yan yield from the much larger number of events 
and expected J/ip/DY ratio. The latter is measured in 
proton-nucleus collisions [28] by NA50 (at 400-450 GeV 
and scaled to 158 GeV) and is corrected for the J/ip "anoma- 
lous suppression" measured in In-In collisions [29]. A 10 % 
relative systematic error is applied to the resulting normal- 
izations, to account for uncertainties in the J/ip anomalous 
suppression and normal absorption patterns. 
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Charm : 2.77±0.27 
Fit x 2 /NDF: 1.1 
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Fig. 13. Signal dimuon mass distribution measured from the 
4000 A (top) and 6500 A (bottom) data samples, compared 
to the superposition of Drell-Yan dimuons (dashed line) and 
muon pairs from open charm decays (dotted line), scaled up 
with respect to the expected yields. 



4.2 Data versus expectations in dimuon mass and 
offset 

Figure 13 compares the signal dimuon mass distributions 
obtained from the 4000 A and 6500 A data samples, inte- 
grated over collision centrality, with the sum of the Drell- 
Yan and open charm contributions. With respect to the 
expected normalizations, described in the previous para- 
graphs, these contributions must be scaled up (by the val- 
ues quoted in the figure) so as to provide the best descrip- 
tion of the measured signal spectrum in the dimuon mass 
window 1.16 < M < 2.56 GeV/c 2 . Within errors, the two 
data samples give perfectly compatible results. A global 
fit gives scaling factors of 1.26 ± 0.09 for Drell-Yan and 
2.61 ± 0.20 for open charm, with x 2 /ndf = 1.02. Further- 
more, essentially the same numerical values are obtained if 
the analysis is redone only selecting events with a match- 
ing x 2 below 1.5 (instead of 3). As previously observed 
by NA38 [2] and NA50 [4], a significant excess of IMR 
muon pairs is observed, which can be well accounted for 
by increasing the charm normalization. 

The big advantage of NA60, with respect to the dimuon 
measurements made by all other heavy-ion experiments, 
is the availability of the dimuon weighted offset variable, 
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Prompt: 2.39±0.13 
Charm : 1 .23±0.25 
Fit x 2 /NDF: 0.5 
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Charm : 1.06±0.22 
Fit x 2 'NDF: 1.5 
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Fig. 14. Same as previous figure but for the dimuon 
weighted offset distributions of the dimuons in the mass range 
1.16-2.56 GeV/c 2 . 



which provides complementary information ideally suited 
to distinguish prompt dimuons from muon pairs stemming 
from displaced decay vertices. Figure 14 shows the dimuon 
weighted offset distribution for the signal dimuons in the 
mass range 1.16-2.56 GeV/c 2 , for the 4000 A and 6500 A 
data samples, compared to the sum of the two contribu- 
tions: prompt dimuons and open charm decays, scaled to 
provide the best fit to data. 

The shape of the prompt dimuon distribution was built 
using the measured dimuons in the J/ip and (j) peaks, 
where the non-prompt signal contributions are less than 
1%. The shape of the open charm distribution was de- 
fined using the muon pairs from the overlay Monte Carlo 
simulation, including the additional smearing needed to 
reproduce the measured 3/i/j and <fi distributions (see [16], 
section 8.4.1 for details). 

The excess dimuons are clearly concentrated in the 
region of small dimuon offsets, excluding the possibility 
that they are due to open charm decays. The best de- 
scription of the measured distribution is obtained when 
the prompts contribution is scaled up by more than a fac- 
tor of two with respect to the expected Drell-Yan yield, 
while the open charm contribution is compatible with the 
yield assumed by NA50 to reproduce the IMR spectra in 
p-A collisions [4]). Fig. 14 underlines the fact that the 
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dimuon weighted offset is a really efficient discriminator 
between prompt and charm dimuons. On the other hand, 
due to the extreme similarity of the mass distributions of 
the excess and charm dimuons (see Fig. 16), it is clear 
that mass alone is completely unable to separate them 
and, therefore, perfectly accommodates the assumption 
that the excess is due to open charm, from the purely 
statistical point of view (see Fig. 13). 

A global fit of both data sets provides scaling factors 
of 2.29 ± 0.08 and 1.16 ± 0.16, for the prompt and open 
charm contributions, respectively (with % 2 /ndf = 1.52). If 
the analysis is redone only using events with matching \ 2 
below 1.5, instead of 3, the results remain the same within 
the statistical errors. If the open charm yield is restricted 
to be within 10% of the nominal value rather than left 
free in the fit, the scaling factors become 2.43 ± 0.09 and 
1.10 ± 0.10, respectively, i.e. again the same within the 
statistical errors. 

Since the tail of the offset distribution lacks the statis- 
tics needed to define the open charm normalization dif- 
ferentially in bins of mass, px & n d centrality, its scale is 
kept fixed to the value determined from the integrated 
sample, corresponding to a full phase space cross section 
<j c5 = 9.5 /ib/nucleon with 14% statistical and 15% sys- 
tematic errors accounting for the uncertainty of the ex- 
pected Drell-Yan contribution. This is the value needed 
to describe the measured data and does not depend on 
the assumed reference cross section. However, the extrap- 
olation of the cross section to full phase space depends 
on the kinematical distributions of the simulated c and c. 
The associated uncertainty is not accounted in the quoted 
systematic error. For instance, we would obtain a value 
19% smaller if we would use the conditions of the NA50 
analysis: Pythia 5.7, different parton densities and with- 
out nuclear modifications, effective cc — > branching 
ratio of 0.97% instead of 0.84%, etc. 



4.3 Kinematic properties of the excess dimuons 

The excess dimuons are defined as the statistical difference 
between the total yield and the sum of fitted charm and 
nominal Drell-Yan. Unbiased physics insight requires cor- 
recting the measured excess for reconstruction efficiencies 
and detection acceptances. The acceptances were calcu- 
lated by Monte Carlo simulations, in dimuon mass and px 
bins, for each of the two data samples, assuming a uniform 
cos #cs distribution and the same rapidity distribution as 
that of the Drell-Yan dimuons (approximately Gaussian 
with a ~ 1). Assuming a Gaussian rapidity distribution 
of width 1.5 gives 20% less excess dimuons, while a 0.5 
width leads to 40 % more excess dimuons. After check- 
ing that the 4000 A and 6500 A data sets give statisti- 
cally compatible results, they were analyzed together. In 
this way we obtain a two-dimensional distribution of the 
acceptance-corrected excess as a function of mass and px- 
Figure 15 shows the ratio between the excess dimuon 
yield (corrected for acceptance with the assumptions just 
mentioned) and the expected Drell-Yan yield (directly tak- 
en from the generator), in the mass range 1.16-2.56 GeV/ c 2 
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Fig. 15. Ratio between the excess and the Drell-Yan dimuon 
yields, after acceptance correction, versus centrality (open cir- 
cles). The excess per squared number of participants is also 
shown, in arbitrary units (filled circles). 

as a function of Y part . The smaller error bars represent the 
statistical errors while the larger ones represent the sum, 
in quadrature, of statistical and systematic uncertainties 
of the fitted Drell-Yan and open charm normalizations. We 
observe that the yield of excess dimuons, per Drell-Yan 
dimuon, increases from peripheral to central In-In colli- 
sions. Furthermore, we see that even the most peripheral 
event sample has a considerable yield of excess dimuons 
(essentially identical to the yield of expected Drell-Yan 
dimuons). 

To gain further insight into the properties of these ex- 
cess dimuons, we have calculated the ratio between their 
yield and Y 2 art . The excess /DY was scaled by the ratio be- 
tween the number of binary nucleon-nucleon collisions and 
the squared number of participant nucleons, N co \\/Np alt , 
calculated for each In-In centrality bin using the Glauber 
model. The resulting excess / N^ &y:t ratio, also plotted in 
Fig. 15, shows a slight decrease with N part . This implies 
that the increase of the excess with centrality is some- 
where in between of linear and quadratic in Y part . 

The mass distribution of the excess dimuons, after ap- 
plying the corrections for acceptance and reconstruction 
efficiencies, is shown in Fig. 16. The mass spectra of open 
charm and Drell-Yan pairs are plotted for comparison. 
The shapes of these spectra are directly taken from the 
generator. The yields reflect the fit values discussed be- 
fore; the bands correspond to the associated systematic 
errors. Clearly, the mass spectrum of the excess drops off 
much more steeply with mass than that of the Drell-Yan 
pairs. In contrast, the slopes are nearly the same for the 
excess and open charm. This explains why the excess, al- 
ready seen by the NA38/NA50 experiment, was found to 
be quite fairly described by any source, be it prompt or 
delayed, with a similar mass distribution as open charm, 
when using the mass spectrum only as a discriminator [6] . 

Fig. 17 summarizes the different aspects of the trans- 
verse momentum distributions of the excess dimuons. The 
top plot shows the excess /DY ratio as a function of the 
dimuon px, clearly indicating that the process responsible 
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Fig. 16. Mass distribution of all three contributions to the 
IMR spectrum, corrected for acceptances and efficiencies as 
described in the text. From top to bottom at the highest mass 
bin: Drell-Yan (as expected from the measured number of J/ip 
dimuons) , excess (the difference between all prompt and Drell- 
Yan pairs) and open charm (with scaling factor of 1.16). 

for the production of the excess dimuons is significantly 
softer than the Drell-Yan dimuons. While in the lowest pt 
bin there arc around 3.5 times more excess dimuons than 
expected Drell-Yan dimuons, this ratio drops to 0.5 at 
high pt- The px spectra themselves are shown in Fig. 17- 
middlc, in 3 consecutive dimuon mass windows: 1.16-1.4, 
1.4-2.0 and 2.0-2.56 GeV/c 2 . Again, a significant devia- 
tion from the behaviour of Drell-Yan pairs is observed. 
The three pt spectra are all different from each other, 
while the px spectra and the mass spectra of Drell-Yan 
factorize: the primordial fcr=0.8 GeV/c characterizing the 
Gaussian distribution is independent of mass in the region 
measured, i.e. from 3 to >10 GeV/c 2 [30]. 

Finally, Fig. 17-bottom shows the same data in terms 
of dimuon transverse mass (rnr = yjv\ + m2 ) spectra, for 
the same three dimuon mass windows defined before. All 
spectra are essentially exponential. However, a steepening 
is observed at very low ttit in the lowest mass window, 
which finds its counterpart in all my spectra observed in 
the low mass region below 1 GeV/c 2 [31], but seems to 
be switched-off in the upper two mass windows as seen 
here. The phenomenon is outside any systematic errors 
as discussed in [31], but has so far not found a convinc- 
ing physical interpretation. Ignoring the Iow-tot rise, the 
data can be fit with simple exponentials 1/pt dN/dpr = 
1/rriT dN/drriT ~ exp{—rriT /T) over the complete Ar- 
range (lines in Fig. 17- middle and bottom), resulting in the 
following respective T c g values: 189 ±15 (stat)±4 (syst), 
197 ± 13 ± 2, and 166 ± 17 ± 4 MeV. The systematic errors 
are dominated by the uncertainties of the Drell-Yan and 
open charm contributions. If the fit is instead restricted to 
Pt > 0.5 GeV/c, consistent with [31] to exclude the rise at 
low-my, the T c g values slightly (but hardly significantly) 
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Fig. 17. Ratio between the excess and the Drell-Yan dimuon 
yields, after acceptance correction, versus dimuon transverse 
momentum (top). Transverse momentum (middle) and trans- 
verse mass (bottom) spectra of the excess dimuons, in three 
dimuon mass windows: 1.16-1.4, 1.4-2.0 and 2.0-2.56 GeV/c 2 . 



rise to 199 ± 21 ± 3, 193 ± 16 ± 2 and 171 ± 21 ± 3 MeV, 
respectively. 



5 Discussion 

The central results of the present paper are connected to 
two essentially independent physics issues: 

— The observed yield of muon pairs from D meson decays 
leads to <r c5 = 9.5 ± 1.3(stat) ± lA(syst) yb (the sys- 
tematic error does not reflect the uncertainty related 
to the kinematic distribution of the c and c quarks) 
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and is compatible with the charm production cross sec- 
tion deduced from the IMR dimuon data measured by 
NA50 in p-A collisions. Charm production in A-A col- 
lisions is not enhanced relative to expectations, and it 
therefore cannot be made responsible for the dimuon 
enhancement seen by NA38 and NA50. 
— The dimuon enhancement in the IMR as observed be- 
fore in Pb-Pb collisions also exists in In-In collisions. 
It is now experimentally proven to be solely due to 
the prompt component. The dimuon excess has been 
statistically isolated by subtracting the Drell-Yan and 
open charm contributions from the total. It is on about 
the same level as Drell-Yan (and charm) and increases 
significantly from peripheral to central collisions rel- 
ative to Drell-Yan and N part . Both its mass and px 
spectra show a much steeper fall-off than Drell-Yan. 
Moreover, the pt spectra depend on mass and do not 
show the factorization between mass and px charac- 
teristic for Drell-Yan. Conversely, fits to the essentially 
exponential mr-spectra lead to inverse slope parame- 
ters T e flf, of about 190 MeV, which do not depend on 
mass, within the (relatively large) errors. 

The prime contender for the interpretation of the ex- 
cess is thermal radiation. The remainder of this section 
places the results reported in this paper into a wider con- 
text, relating them to other experimental results from 
NA60 and to the latest theoretical predictions on thermal 
radiation in the IMR. 

NA60 has also studied dimuon mass and px spectra in 
the low mass region (LMR) M<1 GeV/c 2 [31,32,33]. The 
total mass spectrum, unifying the results from the LMR 
analysis and from the IMR data in Fig. 16, is shown in 
Fig. 18-top. The LMR results correspond to the integral 
of the pr-differential acceptance-corrected mass spectra 
published previously [32,33]; a cut px > 0.2 GeV/c is ap- 
plied to avoid the very large errors in the region of low 
acceptance [32,33]. The mass spectrum is absolutely nor- 
malized in the LMR region as described in [32]; the IMR 
data from Fig. 16 have independently been normalized fol- 
lowing an analogous procedure. Recent theoretical results 
on thermal radiation from two major groups working in 
the field are included for comparison [34,35], calculated 
absolutely (not normalized relative to the data); a further 
result [36] exists, but is left out here due to the lack of final 
normalization. The general agreement between data and 
model results both as to spectral shapes and to absolute 
yields is most remarkable, supporting the term "thermal" 
used throughout this paper. The strong rise towards low 
masses reflects the Boltzmann factor, i.e. the Planck-like 
radiation associated with a very broad, nearly flat spec- 
tral function. Even this part is well described by [34], due 
to the particularly large contribution from baryonic inter- 
actions to the low-mass tail of the p spectral function in 
this model. Higher up in mass, the p pole remains visible, 
followed by a broad bump in the region of the <fi. This is 
described in [34] as in-medium broadening of a small frac- 
tion of the (f> (caution should, however, be presently taken 
on that, since corrections for the resolution function of the 
NA60 apparatus are still under investigation) . In the IMR 
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Fig. 18. Top: Acceptance-corrected mass spectra of the 
excess dimuons for the combined LMR/IMR. Errors in the 
LMR part are statistical; the systematic errors are mostly 
smaller than that. Errors in the IMR part are total er- 
rors. The theoretical model results are labeled according to 
the authors Hees/Rapp [34] (EoS-B + option is used) and 
Renk/Ruppert [35]. Bottom: Inverse slope parameter T e ff of 
the excess mr-spectra vs. dimuon mass [31,32]. For the LMR 
data M<1 GeV/c 2 (triangles), Drell-Yan is not subtracted 
(would decrease the values only within the error bars [31]). 
The IMR data (closed circles) correspond to the present work. 
Open charm is subtracted throughout. The bands show the in- 
verse slopes for the Drell-Yan and open charm contributions 
as provided by Pythia. 

region above 1 GeV/c 2 , the description is good for both 
scenarios (only available up to 1.5 GeV/c 2 for [34]). 

Both for the LMR and IMR, the mr-spectra are pure 
exponentials (at least for px >0.4 GeV/c), consistent with 
the thermal interpretation [35,36]. Apart from the abso- 
lute scale, they can therefore be described by one single pa- 
rameter, the inverse slope T e g extracted from exponential 
fits to the data. The combined results for T c ff in the LMR 
and those reported in this paper are shown in Fig. 18- 
bottom [31]. For M < 1 GeV/c 2 (triangles), a correction 
for Drell-Yan pairs is not done, due to their small contri- 
bution [31], the intrinsic uncertainties at low masses [34] 
as well as the inability of Pythia to generate Drell-Yan in 
this region. The square points correspond to the exten- 
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sion of the LMR analysis up to M = 1.4 GeV/c 2 which 
did not account for the systematic errors of the Drell-Yan 
and open charm contributions. One should note that the 
square points and circles are not statistically independent, 
since the two analyses were performed on overlapping data 
samples. The inverse slopes for the Drell-Yan and open 
charm contributions are shown for comparison. The dif- 
ference to the excess data is most remarkable. 

Below 1 GeV/c 2 , the inverse slope parameters T e g are 
not at all independent of dimuon mass, but monotonically 
rise with mass from the dimuon threshold, where T e fj is 
~180 MeV, up to the nominal pole of the p meson, where 
T e ff is ^250 MeV. This is followed by a sudden decline to 
the level of T c g ~190 MeV reported here. That decline be- 
comes even more steep, jump-like, if the slope parameters 
T e ff are corrected for the contribution of the freeze-out 
p [33]. The initial rise is consistent with the expectations 
for radial flow of a hadronic source (here tt + tt~ — > p) de- 
caying into lepton pairs. However, extrapolating the lower- 
mass trend to beyond 1 GeV/c 2 , a jump of about 50 MeV 
down to a low-flow situation is extremely hard to reconcile 
with emission sources which continue to be of dominantly 
hadronic origin in this region. Rather, the sudden loss of 
flow is most naturally explained as a transition to a qual- 
itatively different source, implying dominantly early, i.e. 
partonic processes like qq — > p + pT for which flow has not 
yet built up, at least at SPS energies, due to the "soft 
point" in the equation-of-state. This may well represent 
the first direct, i.e. data-based evidence for thermal ra- 
diation of partonic origin, overcoming parton-hadron du- 
ality for the yield description in the mass domain (see 
below). The observed slope parameters T e ff ^190 MeV 
are then perfectly reasonable, with a purely thermal in- 
terpretation without much flow, reflecting the averaging 
in the space-time evolution of the fireball between the ini- 
tial temperature T 4 ~220-250 MeV (at the SPS) and the 
critical temperature T c ^170 MeV. 

Theoretically, the NA50 IMR dimuon enhancement [6] 
was successfully described as thermal radiation based on 
parton-hadron duality, without specifying the individual 
sources [7]. However, the same approach is not any longer 
appropriate for the NA60 data. The extension of the uni- 
fied LMR and IMR results over the complete M-p T plane 
places severe constraints on the dynamical trajectories of 
the fireball evolution, allowing for more detailed insight 
into the origin of the different dilcpton sources on the ba- 
sis of radial flow, sensitive to the time ordering of the 
sources. Indeed, all present scenarios [34,35,36] explicitly 
differentiate between hadronic (mostly 47r) and partonic 
contributions in the IMR. The partonic fraction ranges 
from 0.65 for [34] (option EoS-B+ as used in Fig. 18-top) 
to "dominant" in [35,36]. However, due to remaining un- 
certainties in the equation-of-state, in the fireball evolu- 
tion and in the role of hard processes [34] , a quantitative 
description of the very sensitive inverse slope parameter 
T e ff in Fig. 18-bottom is only slowly emerging. In partic- 
ular, the more recent results from the authors of [35,36], 
while very encouraging, are still preliminary and have not 
yet been formally published in their final form. A system- 



atic comparison of several model results to the data in 
Fig. 18-bottom is therefore presently not possible. 



6 Conclusions 

The dimuon excess in the mass region M>1 GeV/c 2 seen 
in high-energy nuclear collisions before has now been pro- 
ven to be of prompt origin. Its properties, differing from 
those of Drell-Yan pairs in many ways, suggest an inter- 
pretation as thermal radiation. If linked to supplementary 
information on dimuon excess production in the mass re- 
gion <1 GeV/c 2 , all indications favor an early, i.e. a domi- 
nantly partonic emission source. Present theoretical mod- 
elling, though still under development, supports our inter- 
pretation. 
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Appendix A : Combinatorial background 

Our aim is to pick pairs of muons from different events in 
such a way that after applying the trigger conditions they 
reproduce the observed like-sign spectra, both in shape 
and in absolute normalization. The problem arises from 
the fact that the data used as pool for the single muon 
sample is already affected by the trigger conditions which 
induces correlations between the registered muons. We 
will now describe the procedure used to account for these 
correlations and to obtain the "unbiased" single muon 
pools. 

We will denote by P + and P~ the average numbers 
of triggerable muons of positive and negative charge, re- 
spectively, in a single interaction (the numerical values 
of these probabilities are always <C 1 since the probabil- 
ity for pion to produce a triggerable muon is ~ 10~ 3 ). If 
we neglect the correlation induced by charge conservation 
between the numbers of positive and negative hadrons (a 
very reasonable assumption in the case of high multiplicity 
heavy-ion collisions) , the number of muon pairs of differ- 
ent charge combinations observed in N collisions will be 

N++ = NP++ = NP+P+/2 
N = NP— = NP-P-/2 (A-l) 
N+- = NP+- = NP+P~ . 

To account for the rejection of the same-sextant muon 
pairs by the trigger, we decompose P + (and P~) in the 
contributions from the different sextants, P + = ^2 i=1 pj 
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and rewrite Eqs. (A-l) excluding the rejected combina- 
tion: 



p ++ = 



E 



Pt Pi 



6 



p 



+- 



Erf 



= p-p- 



p+p- 



-E* 




E 



pi Pi 



It is worth noting that Eqs. (A-2) imply: 



4P 4 



p— = (p + p7 - P'Pi 



Y^pfpj (pip] 

i¥=3 



ptPi) (A-3) 



The well-known equation A H = 2v N++N is a par- 
ticular case of Eq. (A-3) when its right-hand side vanishes. 
Rewriting the latter as 



E 



Eft + ft 



Eft 



_2 _2 




(A-4) 



we see that it vanishes when the pf /pj ratios are the same 
for all sextants (in more general terms, this ratio should be 
constant over the whole phase space) . This is not the case 
in NA60 because the dipole magnet breaks the symmetry 
of the azimuthal distribution of the produced particles, in 
a charge dependent way. Therefore, in order to evaluate 
the CB in NA60, we need to compute the single muon 
probabilities and explicitly account for the exclusion of 
the same-sextant dimuons. Since the number of same-sign 
pairs with muons in the sextants i and j is, for the positive 



case, N^j 
Pi as 



P+P+/2, with p+ 



NpJ' , we can extract 



(Pi 



(A-5) 



and then average over all possible {j, k} combinations. 
Once the values of pf and p~ are known, we can determine 
the fractions of positive and of negative muons, no longer 
being biased by the trigger condition: 



R+ = 



Epi 



Epi + EPi 



R- = 1 - W 



(A-6) 



We can, then, build the artificial CB spectra according to 
the following procedure: 

i. Select randomly the charges of the two muons, accord- 
ing to the probabilities given by Eqs. (A-6). 



ii. Select randomly the sextant of each muon, according 
to the weights p\ and p~ , restarting from step 1 if the 
two selected sextants happen to be the same. 

iii. Randomly pick two muons, of charges and sextants 
as previously selected, from the single muon samples 
built out of the measured like-sign dimuon events. If 
the two selected muons have more than one match to 
VT tracks (mi and mi matches, say), then we build 
all possible (mi x mi) matched dimuons and apply to 
each of them the selection cuts applied to the measured 
events. 

The normalizations of these artificial CB samples are fixed 
by 



^2 pt pj 



Ar ++( - } 4E 



pi^pf-" 



(A-7) 

so that the generated like-sign dimuon spectra reproduce 
the corresponding measured spectra. 

Special care must be taken in what concerns the off- 
sets of the muons in the "mixed" pairs. In order to repro- 
duce the offset distribution of the measured combinatorial 
muons, we first randomly assign the vertex of one of the 
two events participating in the mixing to be the vertex of 
the generated event. Then, we modify the intercept pa- 
rameters of the muon from the other event so that with 
respect to this vertex it retains the same offset as it had in 
its own event, with respect to its vertex. The accuracy of 
the method schematically depicted in Fig. 19 can be ap- 
preciated in Fig. 8, which compares the dimuon weighted 
offset distributions of the generated and measured like- 
sign muon pairs. 




Event 1 



Mixed event 



Event 2 



Fig. 19. Schematic explanation of the method used to build 
the muon offsets for mixed events. 



Appendix B : Fake matches background 

Our aim is to estimate the probability for a given muon 
from the MS to be wrongly matched with the tracks in the 
VT and to use it in building the artificial fake dimuons 
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sample which reproduces both in spectral shape and in 
normalization the FB contributing to data. 

Let e be the probability that the correct match is 
present in the set of all found matches for a muon with 
given kinematics (regardless on the number of fake matches 
and their matching x 2 's). Notice that the correct match 
may be missing not only due to track reconstruction inef- 
ficiency or the choice of the x 2 cut value, but also because 
the muon was produced in a decay or interaction down- 
stream of the VT, in which case its track is simply absent. 
Let us denote by <fi n (n > 0) the probability for a given 
muon to have n fake matches in events with given char- 
acteristics such as interaction sub-target, multiplicity, etc. 
If we designate v as the average number of fake matches, 
then <f> n is Poisson distributed with mean v. Our deriva- 
tions do not require, however, any assumption about its 
distribution. 

These quantities can be extracted for each muon with 
arbitrary precision using a event mixing technique in the 
following way: one applies the usual matching procedure 
but tries to associate the MS muon of one event with the 
VT tracks of many different events with same character- 
istics (i.e. multiplicity, interaction sub-target, etc.). From 
each such event one gets a set of a priori fake matches 
with the same <fi n and \ 2 distributions as for the fakes in 
the real data. The latter distribution, B(x 2 ) (scaled down 
by the number of tried events per muon), after being sub- 
tracted from the x 2 distribution of the real data, provides 
the spectrum for the correct matches, S^X 2 ) (both B(x 2 ) 
and S(x 2 ) are shown in Fig. 6). Then, the probability for 
a given muon to have n matches (either one correct and 
n — 1 fakes or all n fakes) can be written as: 



P(n) = e<f> n -i + (1 - e)(fi n 



(B-l) 



and the probability that the correct match is present in 
this n-plet is 



P pr (n) = ecj> n -i/P(ri) 



(B-2) 



Provided that the correct match is present in this set 
of n matches, Bayes' theorem states that the fraction of 
times Wb(n\pr) in which it will have the smallest (best) 
X 2 is equal to the ratio of the probability for configura- 
tion {xi corr , ...} to the sum of probabilities for all con- 
figurations with given x 2 values (i.e. {xlj a ke> xl.corr., ■■■}, 
etc.). Expressing the probability of given arrangement of 
X 2 values of n matches as the product of the probabilities 
for each x 2 , we can write: 



W b (n\pr) = 



S(xl)I\UB(x 2 ) 



Ri 



Eti s(x 2 )U^b(x 2 ) E"=i R i 



(B-3) 



with R, = S(x 2 )/B(xf)- From Eqs. (B-l-B-3) we get the 
probability for the best one of n matches being the correct 
one: 



W b (n) = P pr (n)W b (n\pr) 



ei?i 



(B-4) 



It is the presence of the e in Eq. (B-4) which makes it 
difficult to estimate the FB in the best matches spectra 
using the data only: the probability of the correct match 
being present in the set of the matches strongly depends 
on the kinematics of the muon. Even for muons with sim- 
ilar kinematics it differs for those coming from the inter- 
action point and those originating in tt or K decays. 

Consider now a pair of muons, each with its own ti 

(i) 

and <fin , i — 1,2. For the sake of generality consider two 
extreme possibilities: the case in which the probability of 
finding the correct match for the first muon is not corre- 
lated with that of the second muon and the case in which 
the correct matches are present or absent always together, 
with common probability e. In the first case the probabil- 
ity of finding n and k matches respectively, similarly to 
Eq. (B-l), can be written as: 



P(n,k)= e 1< fi n 1 l 1 + (1 - er)^ 1 



e 2 . 



i + (l-e 2 ) 



while in the second, full correlation, case we have 



P(n,k) = e$l 1 $l 1 + (1 



(B-5) 



(B-6) 



Taking into account the identities n<fi n - 



Eo° n<j> r , 



1, the average number of 



matched dimuons for given pair, W — J2 n fe=i n kP(n, k), 
is equal to 



W 



e 2 v\ + v\V2\ 



(B-7) 



for the case of absence of the correlation between the cor- 
rect matches and 



W — e + [e(vi + y 2 ) + ^1^2] 



(B-8) 



for the correlated case. Note that in Eqs. (B-7-B-8) the 
expression in square brackets (involving the average num- 
ber of fake matches per muon, v) gives the average number 
of fake dimuons which is what we want to reproduce by 
the event mixing technique. 

We thus arrive at the following procedure for a given 
pair of muons from the MS (including those which have 
no matches) 

i. For each muon of the pair we generate the "mixed 
fakes" by selecting matches from the same number of 
tracks as in the event where the pair comes from, but 
picking the tracks from other events with similar char- 
acteristics. 

ii. Combine all mixed fakes of the first muon with all orig- 
inal matches (if any) of the second one and vice- versa. 
The probability of obtaining this way n x k (including 
the cases of n or k = 0) a priori fake dimuons is 

00 00 

F(n, k) - p ( l > k ) + ^ E P ( n > 

1=0 1=0 

which, after substitution of Eq. (B-6), leads to 



F(n,k) 



e + (i 



4 2) 



e^r + CL-e)^ 1 ) (B-10) 
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in the uncorrelated scenario and 

F(n,fc) = ek 1 )^ 1 + ^i 2) 



(B-ll) 



for the correlated one. Averaging over all possible {n, k} 
combinations (i.e. performing these two steps many 
times with different events using the same muon pair) 
we get for the average number of "mixed fake dimuons" , 
W\, defined as Y^k=i n kF(n, k), 



W\ = eii>2 + £2^1 + v 2 
Wi = e(vi + v 2 ) + 2viv 2 



(B-12) 
(B-13) 



for the uncorrelated and correlated cases, respectively. 

Note that these numbers reproduce the fake dimuons 
contribution (in [ ]) of Eqs. (B-7) and (B-8), respectively, 
except for an extra factor 2 in the term corresponding 
to both matches being fake. This double counting can be 
easily removed by combining the mixed fake matches of 
the first muon with those of the second one and counting 
these dimuons with a negative sign, thus obtaining the 
needed W 2 = —v\v 2 contribution. 

This algorithm does not require any explicit determi- 
nation of e's of <j>nS. For each pair of muons large amounts 
of "mixed" fake dimuons are generated with small weights 
(W\ — W 2 )/N, where N is the number of different events 
matched to the same muon pair, thus smoothing the bin 
to bin fluctuations. 
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